x = gaminv((0.005:0.01:0.995),200,20);
y = gampdf(x,200,20);
y1 = normpdf(x,2000,200);
plot(x,y,'-',x,y1,'-.')
